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ABSTRACT 

We have used the sub-miUimeter array to image the molecular envelope 
around IRC+10420. Our observations reveal a large and clumpy expanding enve- 
lope around the star. The molecular envelope shows a clear asymmetry in ^^CO 
J=2-l emission in the South- West direction. The elongation of the envelope is 
found even more pronounced in the emission of ^^CO J=2-l and SO Jk=65-54. 
A small positional velocity gradient across velocity channels is seen in these lines, 
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suggesting the presence of a weak bipolar outflow in the envelope of IRC+10420. 
In the higher resolution ^^CO J=2 1 map, we find that the envelope has two 
components: (1) an inner shell (shell I) located between radius of about l"-2"; 
(2) an outer shell (shell II) located between 3" to 6" in radius. These shells repre- 
sent two previous mass-loss episodes from IRC+10420. We attempt to derive in 
self-consistent manner the physical conditions inside the envelope by modelling 
the dust properties, and the heating and cooling of molecular gas. We estimate 
a mass loss rate of ^9 10~^ Mq yr~^ for shell I and 7 10~^ Mq yr~^ for shell II. 
The gas temperature is found to be unusually high in IRC+10420 in comparison 
with other oxygen-rich envelopes. The elevated gas temperature is mainly due to 
higher heating rate, which results from the large luminosity of the central star. 
We also derive an isotopic ratio ^^C/^^C = 6. 

Subject headings: circumstellar matter: — ISM: molecules — stars: AGB and 
post-AGB — stars: individual (IRC+10420) — stars: mass loss 

1. Introduction 

Yellow hypergiants are rare objects which are post-red supergiants rapidly evolving in 
blue- ward loops in the Hertzsprung- Russell diagram (de Jager 1998). A few of such ob- 
jects are found to have prodigous mass loss, leading to the formation of thick circumstellar 
envelopes. The study of the circumstellar envelope around them is very important to pro- 
vide insight into the mass loss process from these massive stars and the evolution of stars 
themselves. 

The star IRC+10420 is classified as yellow hypergiant (de Jager 1998). Based on spec- 
troscopic monitoring, the star initially classified as spectral type F8 I^ in 1973 (Humphreys 
et al. 1973), has evolved into spectral type A in 1990 (Oudmaijer et al. 1996, Klochkova et 
al. 1997). Strong emission lines such as the Balmer series and Ca II triplets are seen in its 
optical spectrum (Oudmaijer 1998). Broad wings seen in Ha and Ca II triplets suggest a 
large outfiow velocity in a massive wind close to the stellar surface. High resolution optical 
images of IRC+10420 (Humphreys et al. 1997) revealed a complex circumstellar environment 
with many features such as knots, arcs or loops. Lacking kinematic information, however, 
the real three dimensional structure of the envelope can not be inferred. Recently, through 
integral-field spectroscopy of Ha emission line, Davies et al. (2007) show strong evidence for 
axi-symmetry in the circumstellar envelope. 

The molecular envelope around IRC+10420 is very massive and is a rich source of 
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molecular lines (Quintana-Lacaci et al. 2007). The mass loss rate estimated from CO 
observations is at the upper end seen toward evolved stars, up to a few times 10~^ Mq yr~^ 
(Knapp & Morris 1985, Oudmaijer et al. 1996). High denisty tra cers such as ammonia 



(NH3) and HCN lines have been detected ( IMenten fc Alcolea Ill995l ). Molecules typical of 



oxygen-rich envelope such as SO and SO2 (Sahai & Wannier 1992, Omont et al. 1993) 
are known to be present. IRC+10420 is also a strong OH maser source, exhibiting both 
the main lines (1665 MHZ and 1667 MHz) maser and the satelite line (1612 MHz) maser. 
By comparison with typical OH/IR stars, IRC+10420 has the warmest envelope harboring 
O H masers (Humphreys et al . 1997). High angular resolution imaging of the OH masers 



by iNedoluha fc Bowers I ( 119921 ) shows that the maser emissions are distributed in an oblate 
expanding shell measuring V\5 in radius. Recently, Castro-Carrizo et al. (2007) presented 
high angular resolution (^3^^ and l'^4) mapping of the envelope around IRC+10420 in ^^CO 
J=l-0 and J=2-l lines. Their observations reveal an expanding but complex molecular 
envelope with an inner cavity around the central star, implying a recent dramatic decrease 
in the mass loss from IRC+10420. In addition, the shell shows some elongation in the North- 
East to South- West direction, suggesting a non-spherical symmetry in the mass loss process. 
From excitation analysis of -^^CO lines, Castro-Carrizo et al. (2007) show the presence of 
two separate shells with mass loss in the range 1.2 10"^ to 3 1O"^M0 yr-\ Interestingly, the 
gas temperature in shells derived from their model fitting is unusually high, 100 K even at 
the large radial distance of ^2 10^'' cm. 

The distance to IRC+10420 is still uncertain. From analysis of photometric and po- 
larimetric data, Jones et al. (1993) conclude that IRC+10420 is located at a large distance, 
between 4 to 6 kpc. Here we adopt a distance to IRC+10420 of 5 kpc, similar to that adopted 
in Castro-Carrizo et al. (2007). 

In this paper we present high angular resolution observations of the molecular envelope 
around IRC+10420. We also perform detailed modelling to derive in a self-consistent manner 
the dust properties and molecular gas temperature profile in the envelope. 



2. Observations 

We have observed IRC+10420 in the 1.3mm band with the SubMillimeter Array (SMA) 
at two different epochs and in two different configurations as summarized in Table 2. The first 
observation was carried out on 2004 June 2P^, with 8 antennas in an extended configuration, 
giving projected baselines from 27 to 226 m. The weather was good during the observation 
with zenith atmospheric opacity of r225GHz ^ 0.2 at 225 GHz and the system temperatures 
ranged between 200 and 600 K. The second observation was carried on 2005 July 2^^, with 
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Table 1: Basic data of IRC+ 10420 and its molecular envelope 



R.A. (J2000) 19^26^48^03 (1) 

Dec. (J2000) ir2n6'.7 (1) 

Distance 5 kpc (2) 

Angular scale V - 5000 AU - 7.5 10^^ cm 

Luminosity L* '--^ 6 10^ Lq (2) 

Spectral Type mid-A type (3,4) 

Effective Temperature Teff 7000 K (3) 

Stellar radius 1.66 Au 

Systemic velocity Vlsr = 74 kms"-^ (3) 

Expansion velocity Vgxp = 40 kms"-^ (3) 

Mass loss rate M - 6.5 10"^ yr-^ (5) 



References. — (1) SIMBAD database; (2) Jones et al. (1993); (3) Oudmaijer et al. (1996); (4) Klochkova 
et al. 1997; (5) Knapp & Morris (1985) 

7 antennas in a compact configuration, giving projected baselines in the range from 11 to 
69 m. The weather was very good with r225GHz ^0.1 and the system temperatures ranged 
between 100 and 200 K. 

In each observation, the correlator of the SMA was configured in the normal mode with 
two 2-GHz bandwidth windows separated by 10 GHz. This allowed us to cover simulta- 
neously the ^^CO J=2-l line in the upper sideband (USB) the ^^CO J=2-l and SO 65-54 
lines in the lower sideband (LSB). The bandpass of individual antennas was calibrated using 
the quasar 0423-013 for the observation taken in 2004 and Uranus for observation taken in 
2005. To improve the S/N of the data, we smoothed data from the instrumental spectral 
resolution of 0.8125 MHz to 3.25 MHz {i.e. - 4 km s"^). For comparison, the observed lines 
have velocity widths of (~ 80 km s~^). The nearby quasars 1751+096 and 1925+211 were 
observed at regular intervals to correct for gain variations of the antennas due to atmospheric 
fluctuations. The data reduction and calibration were done under MIR/IDl£|. After the gain 
and bandpass calibrations, the visibilities data from both tracks were combined together. 
Imaging and deconvolution were performed under GILDA^ package. 



^http:/ /cfa- www. harvard.edu/~cqi/mircook. html 
^http: / /www.iram.fr /IRAMFR/GILDAS / 
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The 1.3 mm (225 GHz) continuum emission map was derived by averaging the line- free 
channels from both LSB and USB, resulting in a total bandwidth of 3.4 GHz. For the line 
channel maps we subtract the continuum visibilities of the relevant sideband from the line 
emission visibilites. In Table [2] we summarize the observations of IRC+10420 carried out by 
the SMA. 

For the sake of presentation clearity, the center of all the maps were shifted in the 
visibility domain so that the phase reference matches the position of the continuum peak 
emission, i.e. R.A. = 19^26^48^09 a nd Dec. = lV2VlQ".7 b , whic h also corresponds to the 



center of the SiO shell detected by ICastro-Carrizo et al. I (120011 ) . We note that the lack 



of interferometric baselines shorter than 11 m makes our observations blind to extended 
structures of size > 2h" . In order to estimate the amount of flux recovered by the SMA, we 
convolved the channel maps of ^^CO J=2-l and ^^CO J=2-l lines to the spatial resolution 
of 19'^7 and 12^', respectively. We show in Figure 2 the spectra taken from the these maps 
and the single dish observations taken by Kemper et al. (2003) and Bujarrabal et al. (2001) 
at the same angular resolutions. We estimate that the SMA recovered more than 60% of the 
flux in ^^CO J=2 1 line and almost aU the flux in ^^CO J=2~l line. 

Table 2: SMA line observations towards IRC+10420. 



Line 


Rest Freq 


Config. 


Clean beam 


P.A. 


AV 


Icr rms 




(GHz) 




(FWHM, arcsec) 


(deg) 


(km s~"^) 


(mJy/beam) 


i2CO(2-l) 


230.538 


Comp 


3.44 X 3.12 


46° 


4.2 


60 






Comp+Ext 


1.44 X 1.12 


46° 


4.2 


80 


i3CO(2-l) 


220.399 


Comp 


3.52 X 3.21 


63° 


4.4 


50 


30(65-54) 


219.949 


Comp 


3.51 X 3.19 


63° 


4.4 


50 



Results 



3.1. 1.3mm continuum emission 



As shown in Figure 1, the 1.3mm continuum emission from IRC+10420 was well detected 
in our SMA data. The emi ssion is not resolv e d and has a flux density ^ 45 mJy. A previous 
single-dish observation by IWalmsley et al. I ( 119911 ) also at 1.3mm gives a flux of 101 mJy. 
That suggests that the dust emission is extended, with about 55% of the continuum flux 
resolved out by the interferometer. Assuming the central star to have a black body spectrum 
with a temperature of Teff = 7000 K, it contributes only a negligible amount of the flux at 
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1.3 mm of S* ^ 0.1 mJy. Thus most of the 1.3mm continuum emission is produced by the 
warm dust in the envelope around IRC+10420. 

3.2. ^^CO 3=2-1 emission 

In Figure. 3 we show the channel maps of the ^^CO J— 2—1 emission obtained with 
the compact configuration of the SMA. At our angular resolution of ^3'', the envelope 
is well resolved. At both blueshifted and redshifted velocities the emission appears to be 
circularly symmetric, as is expected for a spherically expanding envelope. More interestingly, 
at velocities between 59 - 89 kms~^, centered around the systemic velocity, the envelope is 
clearly elongated in approximately the North-East to South- West direction. By fitting a two- 
dimensional Gaussian to the intensity distribution of the ^^CO J=2— 1 emission at systemic 
velocity Vlsr = 74 kms~^, we obtain a position angle for the ^^CO J=2— 1 envelope of 
PA=70°. The orientation of the CO emission is consistent with the elongation found by 
Castro-Carrizo et al. (2007). The position- velocity (PV) diagrams of the ^^CO J=2-l 
emission along the major (PA ^70°) and minor (PA ^160°) axis of the envelope as presented 
in Figure 7 are typical of an expanding envelope at an expansion velocity of ^40 kms~^. 

We show in Figure. 4 the channel maps of the ^^CO J=2— 1 emission, which is obtained 
by combining the data from both compact and extended configurations of the SMA. The 
^^CO J=2— 1 emission is further resolved into more complex structures. In the velocity 
channels between 59 - 89 kms~^, centered around the systemic velocity, the ^^CO J=2-l 
emission consists of a central prominent hollow shell structure of to 2'' in radius and a 
more clumpy arc or shell located between radii of 3^' to &\ which is more prominent in the 
South West quadrant of the envelope. Close inspection of the velocity channels around the 
systemic velocity Vlsr ^ 74 km s~^ indicates that the central hollow shell-like structure is 
clumpy and shows stronger emission toward the South- Western side. The outer arc or shell 
of emission can be seen more easily in the azimuthal average of the ^^CO J=2— 1 brightness 
temperature as shown in Figure. 13. We can clearly see a central depression within a radius 
of ^ 1" and enhanced emission between the radii of 3'' to 6''. 



3.3. i^CO J=2-l emission 

The channel maps of ^^CO J=2-l emission obtained from the compact configuration 
data at an angular resolution of ^3''. 5 is shown in Figure 5. The emission in the ^^CO J=2-l 
line appears fainter and more compact then seen in the main isotope ^^CO J=2-l line. Like 
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in ^^CO J=2-l emission, near the systemic velocity, the envelope traced by ^^CO J=2-l 
emission is elongated at position angle of PA ^ 70°. The centroid of the emission appears 
to be slightly offset from the stellar position. Also like in ^^CO J=2-l emission, at higher 
velocities in both blue-shifted and red-sfhited parts the ^^CO J=2-l emission again shows a 
roughly circularly symmetric morphology, as would be expected for a spherically expanding 
envelope. The spatial kinematics of the ^^CO J=2-l emission can be more clearly seen in 
the position- velocity diagram (Figure 7) along the major axis (PA ^70°) of the envelope. 
In the PV diagram, there is a small velocity gradient in the velocity range between 70 to 
100 kms~^. Closer inspection of the channel maps suggests that the small velocity gradient 
corresponds to the slight shift of the emission centroid toward the North-East quadrant of 
the envelope in the abovementioned velocity range. 



3.4. SO Jk=65-54 emission 

In Figure. 6 we show channel maps of the SO 65—64 emission obtained from the data of 
the compact configuration. The emission is well resolved in some of the velocity channels. A 
close inspection reveals that the distribution of SO 65—64 emission is different from that of 
^^CO J=2— 1. At extreme redshifted velocities the centroid of the emission is shifted to the 
East, whereas at extreme blueshifted velocities the centroid of the emission is shifted to the 
West. In addition, the SO 65—64 emission in velocity channels centered around the systemic 
velocity is clearly elongated to the South- West. This elongation is very similar to that seen 
in ^^CO J=2— 1 emission at the same angular resolution (see Figure 3). More interestingly, 
at the opposing redshifted velocities of 84 to 102 kms~^(see Figure 6) the emission is very 
compact and shifts to the North-East of the central star. Such positional shift is shown more 
clearly in the position-velocity diagram of the SO Ji^=65-64 emission along the major axis 
(PA ^70°) of the envelope. There is a clearly velocity gradient in the velocity range between 
70 to 90 kms~^, which is similar but more pronounced that that seen in ^^CO J=2-l line. 

Thus, the SO 65—64 emission in particular but also the ^^CO J=2— 1 emission seem 
to better trace the asymmetric structure inside the envelope. The small velocity gradient 
seen in both lines points to the presence of an ejecta in the envelope of IRC+10420. From 
the abovementioned elongation of the envelope and the orientation of the positional velocity 
shift, we estimate that the ejecta is oriented at position angle of PA=70°. We note that 
SO emission is known to be enhanced in AGB and post-AGB envelopes where collimated 
high velocity outflows are present. Examples can be found in the rotten egg nebula OH 
231.8+4.2 (Sanchez Contreras et al. 2000) and VY CMa (MuUer et al. 2007). In both cases 
SO emission is found to be strongly enhanced in the bipolar lobes. The SO enhancement 
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in bipolar outflows could be related to the synthesis of SO in the shocked molecular gas in 
these outflows. Alternatively, the warm environment around supergiants like VY CMa and 
IRC+ 10420 could facilitate the synthesis of SO through chemical pathways such as S + OH 
— > SO + H or HS + O — > SO + H (Willacy & Millar 1997). 

4. Structure of the envelope 

In this section we attempt to build a model for the circumstellar envelope of IRC+10420 
using our newly obtained SMA data and also the large body of observational data on 
IRC+10420, ranging from optical to milimeter wavelengths. Our goal is to better under- 
stand the physical conditions of the molecular gas inside the envelope and also to retrace the 
mass loss history of the central star. The structure of the envelope around IRC+10420 has 
been previously modelled in the work of Castro-Carrizo et al. (2007), by prescribing a radial 
distribution for the temperature in order to derive the gas density and consequently the mass 
loss rate. In our model we take into account the properties of the dust inside envelope and 
the balance of heating and coohng processes. From the discussion in the previous section, 
the envelope can be approximated as spherically symmetric. The physical conditions such as 
kinetic temperature, density in the envelope can then be inferred from detailed calculations 
of the energy balance within the envelope and matching the predictions of the radiative 
transfer model to the observed strength of the CO rotational lines. This procedure repre- 
sents an improvement in comparison to previous models, and allows a better understanding 
of the physical conditions in the envelope of IRC+10420. 

4.1. A simple model for dust continuum emission 

To determine the heating process in the detected molecular shells requires the knowledge 
of the momentum transfer coefficient (or the flux average extinction coeflicient < Q >) 
between the dust particles and the gas molecules. To estimate this quantity, we follow the 
treatment of Oudmaijer et al. (1996) who modelled in detail the spectral energy distribution 
(SED) of IRC+10420. We emphasize here that by following the treatment of Oudmaijer et 
al. (1996), the structure of the molecular envelope as seen in our SMA data and previously 
in the data of Castro-Carrizo et al. (2007) is not considered exphcitly. 

We use the photometric data of IRC+10420 collected by Jones et al. (1993) and Oud- 
maijer et al. (1996), i.e. the 1992 photometric dataset. We also assume that the dust 
consists of silicate particles with a radius of a = 0.05 /im. The dust opacity as a function 



-9- 



of wavelength is taken from the work of Volk & Kwok (1988). Following Oudmaijer et al. 

(1996), we scale the opacity law to match the absolute value of the opacity at 60/im of /^eo/im 
= 150 cm^ g~^. We also a use a specific density for dust particles of p = 2 g cm~^ and 
assume a dust to gas ratio ^ = 5 10~^ as used by Oudmaijer et al. (1996). 

To fit the SED of IRC+10420, especially in the mid-IR region between 5 to 20 /im, 
Oudmaijer et al. (1996) suggest that two dust shells are needed: a hot inner shell with a low 
mass loss rate and an outer cooler shell with a higher mass loss rate. We note that Blocker 
et al. (1999) also reached same conclusion in their attempt to fit the SED of IRC+10420. In 
our model we also use two dust shells: a hot inner shell with a low mass loss rate of M = 8 
10"^ Mo/yr located between 2.6 10^^ cm (70 stellar radh) to 1.8 10^^ cm (486 stellar radh) 
and an outer cooler shell with a higher mass loss rate of M = 1.2 10~^ Mq/jt extending to a 
radial distance of 10^^ cm where the dust particle number density drops to very low values. 
We also infer a stellar luminosity of 6 10^ Lq. The size of the dust shells, the mass loss rates 
and the stellar luminosity are very similar to that derived by Oudmaijer et al. (1996) scaled 
to the adopted distance of 5 kpc. The results of our radiative transfer model are presented 
in Figure 8 and 9. The fiuxes from the model are corrected for an interstellar extinction of 
Av=5 as suggested by Oudmaijer et al. (1996). As can be seen in Figure 8, the results of 
our model are in close agreement with that shown in Oudmaijer et al. (1996) and match 
reasonably the observed SED of IRC+10420. As shown in Figure. 9, the dust temperature 
decreases monotonically from ^1000 K at the inner radius to ^50 K at the outer radius of 
the dust envelope, even though there is a jump in mass loss rate between the inner hot dust 
shell and the outer cooler dust shell at a radial distance of 1.8 10^^ cm. 

From the flux Fx of the radiation fleld at each point in the envelope we can deflne 
the flux average extinction coefficient < Q >, which characterizes the transfer of angular 
momentum from radiation photons to dust particles during the absorption and reemission 
processes of the continuum radiation. Because in our case we do not consider the scattering 
of radiation, the extinction coefl&cient Qx is simply the opacity where A is the wavelength 
of the radiation. 

The fiux average extinction coefficient < Q > changes significantly in the inner part of the 
dust envelope but reaches a constant value in the outer part. Because the molecular shells 
are located in the outer part of the dust envelope, at radii larger than a few times 10^^ cm, 
for simphcity, we will adopt a constant value < Q >=0.025 (see Figure 9). 
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4.2. 



Excitation of CO molecule 



The temperature of the molecular gas at any point in the envelope is determined by 
the balance between cooling due to adiabatic expansion and molecular emissions, and the 
heating, which is mainly due to collisions of molecules with dust grains (Goldreich & Scoville 
1976). As the dust grains stream through the gas under radiation pressure, the grains transfer 
angular momentum to the gas by collisions. The heating term is therefore directly related 
to the drift velocity: 



where rid is the number density of dust particles. The drift velocity is determined by the 
relation: 



and c is the speed of light, < Q > is the flux average extinction coeflicient of dust particles 
estimated in previous section. The same dust to gas ratio of * = 5 10~^ is used to calculate 
the heating rate. 

Generally, in oxygen-rich envelopes the main coolants are H2O and CO molecules. The 
contribution of H2O molecules to the cooling process in IRC+10420 is quite uncertain. We 
note that far-IR emission lines of H2O are not detected in the ISO data (Molster et al. 
2002). In addition, OH masers are seen in a shell of radius 1 to 1.5 arcsec (Nedoluha & 
Bowers 1992). Because OH radical is the photodissociation product of H2O, it is likely that 
H2O molecules can only exist at even smaller radii, i.e. inside the central cavity disccused 
previously. Therefore, we do not expect the H2O molecules to be present in significant 
quantity in the part of the molecular envelope traced in CO emission As a result, we do not 
include the cooling due to the H2O molecule in our model and calculate the cooling solely 
due to CO molecule and its isotopomer ^^CO. 

We use the thermal balance and radiative transfer model published by Dinh-V-Trung & 
Nguyen-Q-Rieu (2000) to derive the temperature profile of the molecular gas and to calculate 
the excitation of CO molecules and predict the strength of CO rotational transitions. We take 
into account all rotational levels of CO molecules up to J=20. The coUisional cross sections 
with H2 are taken from Flower (2001), assuming an ortho to para ratio of 3 for molecular 
hydrogen. The local linewidth is determined by the turbulence velocity and the local thermal 
linewidth. Because the expansion velocity of the gas in the envelope of IRC+10420 is large 
and the observed CO line profiles have smooth slopes at both blueshifted and redshifted 
edges (see Figure 11), we need to assume a constant turbulence velocity of 3 kms~^. We 



(2) 




(3) 



where is the stellar luminosity, Vexp is the expansion velocity, M is the mass loss rate. 
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note that this is higher than normally used (about 1 kms~^) for the envelope around AGB 
stars where the expansion velocity is lower (Justtanont et al. 1994). In our model, the 
spherical envelope is covered with a grid of 90 radial points. For each rotational transition, 
the radiative transfer equation is integrated accurately through the envelope to determine 
the average radiation field at each radial mesh point. We use the approximate lambda 
operator method together with Ng-acceleration to update the populations of CO molecules 
of different energy levels. Once convergence is achieved, we convolve the predicted radial 
distribution of the line intensity with a Gaussian of specified FWHM to produce simulated 
data for comparison with observations. 

As discussed in the previous section, there is a clear evidence for the presence of two 
molecular shells from our SMA data. The size and location of each molecular shell are 
determined from fits to the interferometric data. In our model we adjust both the mass loss 
rates and the relative abundance of CO molecules in the two molecular shells to match both 
the single dish and interferometric data. We find that a relatively low mass loss rate similar 
to that derived by Castro-Carrizo et al. (2007) would result in a too high gas temperature in 
the envelope as the heating term becomes more dominant in comparison to the cooling term. 
We also find that a high abundance of CO of a few times 10~^ as usually used for oxygen rich 
envelopes (Kemper et al. 2003, Castro-Carrizo et al. 2007) would result in CO lines that are 
too strong. Instead, a relative abundance of 10~^ for CO molecules with respect to H2 for 
the CO molecules provides reasonable fit to all the available data. The radial profile of the 
gas temperature is shown in Figure 10. The sudden jump in gas temperature at the inner 
radius of the second shell (^1.85 10^ cm) is due to the change in the mass loss rate. We can 
see that the gas temperature in the envelope of IRC+10420, although lower than assumed in 
Castro-Carrizo et al. (2007), is higher than that typically found in circumstellar envelopes. 
Even at the outermost radius of the envelope, the gas temperature is still around 50 K. 
For comparison, in the envelope around OH/IR stars, the gas temperature at similar radial 
distances from the central star is predicted to be much lower, 10 K or even less (Justtanont 
et al. 1994, Groenewegen 1994). The main reason for the elevated gas temperature in the 
envelope of IRC+10420 is the strong heating due to the enormous luminosity of the central 
yellow hypergiant. The cooling due to molecular emissions is also affected by the lower 
abundance of CO, that we find neccessary to match the data. A summary of the parameters 
used in our model is presented in Table 3. 

As shown in Figures 11-14, the results of our model are in reasonable agreement with 
both the strengths of CO lines measured by single dish telescopes and the radial distribution 
of brightness temperature of J=l and J=2-l lines observed with interferometers. We note 
that except for the extra emission in the blue part of the line profile between 65 and 75 
kms~^, the shape of the J=2-l and J=3-2 is in close agreement with observation. Higher 
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transitions J=4-3 and 6-5 show some evidence of narrower linewidth than predicted by the 
model. One possible explanation is that the expansion velocity in the inner region traced by 
these high lying transitions has slightly lower expansion velocity than in the outer part of 
the envelope. 

By comparing the model predictions with the strengths of ^^CO J=l-0 and J=2-l 
lines observed with the IRAM 30m telescope (Bujarrabal et al. 2001), we derive a relative 
abundance for ^^C0/H2 of 1.5 10~^, or an isotopic ratio ^^C/^^C ^ 6. This ratio is very 
similar to the values ^^C/^^C = 6 and 12, respectively, found in the red supergiants a Ori and 
a Sco (Harris & Lambert 1984, Hinkle et al. 1976). Such a low isotopic ratio suggests that 
IRC+10420 has experienced significant mixing of H burning products to its surface prior to 
the ejection of the material in the envelope. 

The difference between the mass-loss rates and gas temperature profile derived from 
our modelling and that obtained by Castro-Carrizo et al. (2007) can be understood as the 
consequence of the higher gas temperature and higher CO abundance adopted in their work. 
For the same strength of the CO rotational lines, the higher CO abundance can be almost 
exactly compensated by a corresponding decrease in the mass loss rate. 

Our results indicate the presence of two separate shells (shell I and shell II, see Table 
3) with slightly different mass loss rates. With an expansion velocity of 38 kms~^, the time 
interval between the two shells is ^ 200 yrs. The cavity inside shell I also implies a dramatic 
decrease in mass loss from IRC+10420 over the last ^300 yrs. Thus IRC+10420 loses mass 
in intense bursts, separated by relatively quiet periods of a few hundreds years in duration. 



Table 3: Parameters of the model for IRC+10420 envelope. 



Parameters 


Shell I 


Shell II 


Inner radius 


3.5 10^*5 cm 


1.85 10^^ cm 


Outer radius 


1.5 10^^ cm 


5 10^^ cm 


Mass loss rate M 


9 10-^ M0 yr-i 


7 10-4 M0 yr-i 


Expansion velocity Vexp 


38 kms"^ 


38 kms"^ 


Dust to gas ratio ^ 


5 10"=^ 


5 10"=^ 


Momentum transfer coefficient < Q > 


0.025 


0.025 


12CO/H2 


10-^ 


10-4 


13CO/H2 


1.5 10-5 


1.5 10-5 


Turbulent velocity 


3 kms"^ 


3 kms-^ 
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4.3. Caveats 

In our model we consider separately the dust and gas component in the envelope of 
IRC+10420. The dust continuum emission model is used as a simple way to estimate the 
< Q > parameter of dust particles, which is needed for the thermal balance calculation in 
the molecular shells. It turned out that the mass loss rates of the molecular shells between 
9 10~^ M0/yr for shell I and 7 10~^ M0/yr for shell II are lower but quite close to the 
mass loss rate of 1.2 10~^ M0/yr of the cool outer dust shell inferred from the modeUing the 
SED of the envelope. Our approach, although less ad hoc than simply assuming a value of 
the < Q > parameter for dust particles, is not fully self-consistent in the treatment of the 
dust and gas component. It would be desirable in the future to treat both the dust and gas 
component together when high angular resolution data on the dust continuum emission and 
higher J lines of CO become available. 

Another caveat of our model is the assumption that the heating of molecular gas is 
mainly due to collision between dust particles and gas molecules. The high luminosity of the 
central star of IRC+10420 suggests that radiation pressure on dust particles is a possible 
mechanism for driving the wind as assumed in our mode. However, given the high expansion 
velocity and the presence of many small scale features in high angular resolution optical 
images of the envelope (Humphreys et al. 1997), which has been interpreted as jets or 
condensations, other mechanisms such as shocks produced by interaction between higher 
velocity ejecta and the envelope should be considered. The heating due shocks within the 
envelope might contribute to the thermal balance of the molecular gas in IRC+10420. 

5. Comparison with optical imaging data 

The circumstellar envelope around IRC+10420 is known to be very complex, with a 
number of peculiar structures (Humphreys et al. 1997). In the outer part of the envelope 
at radial distance of ^3'' to 6'^ several arc-like features can be seen. These features have 
been interpreted as representing different mass loss episodes of IRC+10420. In addition, 
condensations or blobs are seen closer (about V to 2'' in radius) to the central star. That 
might represent a more recent mass loss episode. These features correspond spatially to the 
two molecular shells (shell I and II) identified in our observations. In addition, Humphreys 
et al. (1997) also identified several broad fan-shaped features to the South- West side of the 
central star, between radius of 0" .5 to 2" . From their higher surface brightness in comparison 
with other parts of the envelope, Humphrey et al. (1997) suggest that these features are 
ejecta moving obliquely toward the observer. As discussed in previous section, the channel 
maps of ^^CO J=2-l and SO Jk=65-54 (Figures 5 and 6) reveal enhanced emissions in 
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blueshifted velocity channels together with a positional shift of the emission centroid also 
to the South- West quadrant of the envelope. Therefore, the spatial kinematics obtained 
from our observations and the properties of these fan-shaped features are consistent with 
the presence of an ejecta in the South- West side of the envelope around IRC+10420. 

6. Summary 

We have used the sub-millimeter array to image and study the structure of the molecular 
envelope around IRC+10420. Our observations reveal a large expanding envelope with a 
clumpy and very complex structure. The envelope shows clear asymmetry in ^^CO J=2- 
1 emission in the South- West direction at position angle PA^70°. The elongation of the 
envelope is found even more pronounced in the emission of ^^CO J=2 1 and SO Jk=65-54. 
A positional shift with velocity is seen in the above emission lines, suggesting the presence 
of a weak bipolar outflow in the envelope of IRC+10420. 

In the higher resolution data of ^^CO J=2-l, we find that the envelope has two com- 
ponents: (1) a inner shell (shell I) located between radius of about l"-2"; (2) an outer shell 
(shell II) between radius 3" to 6". These shells represernt two previous mass loss episodes 
from IRC+10420. 

We constrain the physical conditions inside the envelope by modelling the dust proper- 
ties, the heating and cooling of molecular gas. From comparison with observations we derive 
the size and the location of each molecular shell. We estimate a mass loss rate of ^ 9 10~^ 
M0 yr~^ for shell I and 7 10~^ Mq yr~^ for shell II. The gas temperature is found to be an 
usually high in IRC+10420 in comparison with other oxygen rich envelopes, mainly due to 
the large heating induced by the large luminosity of the central star. 

We also derive a low isotopic ratio ^^C/^^C = 6 for IRC+10420, which suggests a strong 
mixing of processed material from stellar interior to the surface of the star. 

We are grateful to the Sub-millimeter array staff for their help in carrying out the 
observations. The SMA is a collaborative project between the Smithsonian Astrophysical 
Observatory and Academia Sinica Institute of Astronomy and Astrophysics of Taiwan. We 
thank F. Kemper and D. Teyssier for providing single dish CO spectra of IRC+10420. This 
research has made use of NASA's Astrophysics Data System Bibliographic Services and the 
SIMBAD database, operated at CDS, Strasbourg, France. 
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IRC+10420 1.3 mm 




19^26"'48'5 48'0 47'5 



Fig. 1. — 1.3 mm continuum emission map of IRC+10420 obtained with the SMA compact 
configuration. The synthesized beam is 3'^52 x 3'^21 (PA = 54°) as indicated by the eUipse 
in the bottom left corner. Contour levels are every 6 mJy beam~^ 3a). 
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0.05 



50 100 
LSR velocity (km/s) 

Fig. 2.— Spectra of ^^CO J=2-l (top), ^^CO J=2-l (middle) and 80(65-54) {bottom) 
transitions towards IRC+10420. The spectra are shown in main beam brightness tempera- 
ture with a telescope beam as mentioned in top right corner of each box. The dotted lines 
show the spectra of ^^CO J=2-l obtained with JCMT (Kemper et al. 2003) and ^^CO 
J=2-l obtained with the IRAM 30m (Bujarrabal et al. 2001) 
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IRC+10420 ^^C0(2-l) 

' ' I ' ' ' ' I ' ' ' ' I ' ' '_L' ' I I ' ' ' ' I ' ' ' ' I ' ' '_L' ' ' I ' ' ' ' I ' ' ' ' I ' ' '_L' ' ' I ' ' ' ' I ' ' ' ' I ' ' '_L' ' ' I 

21 :: 26 :: 30 :: 34 :: 38 




Fig. 3. — Channel maps of the ^^CO J=2-l hne emission obtained from the SMA compact 
configuration alone. The synthesized beam is 3^^44 x 3'M2 (PA = 46°) and is indicated in 
the bottom left corner of the figure. Contour levels are drawn every 0.6 Jy beam~^ or 1.3 
K 10a). The velocity resolution is 4.2 km s~-^. The cross indicates the position of the 
continuum peak emission. 



-20- 




Fig. 4. — Channel maps of the ^^CO J=2-l hne emission obtained by combining data from 
compact and extended configurations. The synthesized beam is l'^44 x V\12 (PA = 46°) 
and is indicated in the bottom left corner of the figure. Contour levels are drawn every 0.24 
Jy beam"-^ or 3.4 K 3a). The velocity resolution is 4.2 km s~-^. The cross indicates the 
position of the continuum peak emission. 
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IRC+10420 C0(2-l) 




R.A. Offset (arcsec.) 

Fig. 5. — Channel maps of the ^^CO J=2-l hne emission obtained from the SMA compact 
configuration alone. The synthesized beam is 3^^52 x 3'^21 (PA = 63°) and is indicated in 
the bottom left corner of the figure. Contour levels are drawn every 0.2 Jy beam~^ or 0.45 
K 4a). The velocity resolution is 4.4 km s~-^. The cross indicates the position of the 
continuum peak emission. 
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IRC+10420 30(65-54) 
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Fig. 6. — Channel maps of the SO Jk=65-54 hne emission obtained from the SMA compact 
configuration alone. The synthesized beam is 3^^51 x 3'^21 (PA = 63°) and is indicated in 
the bottom left corner of the figure. Contour levels are drawn every 0.2 Jy beam~^ or 0.45 
K 4a). The velocity resolution is 4.4 km s~-^. The cross indicates the position of the 
continuum peak emission. 
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Fig. 7. — Position- velocity diagrams along the major axis, defined at PA=70° (left) and 
perpendicularly along the minor axis (right). Contour levels are every 0.24 Jy beam~^ for 
^^CO J=2-l and every 0.2 Jy beam-^ for ^^CO J=2-l and SO 65-54. 
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Fig. 8. — Spectral energy distribution of IRC+10420. The crosses denote observational data 
from Jones et al. (1993) and Oudmaijer et al. (1996). The solid line represents the model 
SED. 
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Fig. 9. — Dust temperature profile (dotted line) and flux averaged extinction coeflicient < 
Q > (solid line) in the envelope of IRC+10420 obtained from matching the model prediction 
of the SED to the observational data. 
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Fig. 10. — The predicted radial profile of the molecular gas temperature in the envelope of 
IRC+10420. 
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Fig. 11. — Comparison between observations of CO transitions from Kemper et al. (2003) 
and Teyssier et al. (2006) and our model predictions. The telescope beams are 22" for J=l-0, 
19".7 for J=2-l, 13".2 for J=3-2, 10".8 for J=4-3 and 10".3 for J=6-5 line, respectively. 
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Fig. 12. — Comparison between the observed ^^CO J=l-0 and J=2-l lines observed by 
Bujarrabal et al. (2001) and our model predictions (thick sohd line). The telescope beams 
are 22'' for J=l-0 and 12'' for J=2-l line, respectively. 
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Fig. 13. — Azimuthal average of the ^^CO J=2-l brightness temperature (dotted Une) ob- 
tained from SMA data at velocity channel Vlsr = 75 kms~^, and the model prediction (solid 
line) . 
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Fig. 14. — Comparison between the radial distribution of surface brightness (dotted hue) of 
CO J=l-0 and J=2-l observed with the Plateau de Bure interferometer (Castro-Carrizo et 
al. 2007) and our model predictions (solid line). The angular resolution is 2". 86 for the CO 
J=l-0 and 1"A for the CO J=2-l, respectively 



